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ABSTRACT 

The statistical properties of magnetic discontinuities in the solar wind are investigated by 
measuring fluctuations in the magnetic field direction, given by the rotation A9 that the magnetic 
field vector undergoes during time interval At. We show that the probability density function for 
rotations, P{A9), can be described by a simple model in which the magnetic field vector rotates 
with a relative increment AB/B that is lognormally distributed. We find that the probability 
density function of increments, P{AB / B), has a remarkable scaling property: the normalized 
variable x = {AB / B) ■ {At / Ato)^°' has a universal lognormal distribution for all time intervals At. 
We then compare measurements from the solar wind with those from direct numerical simulations 
of magnetohydrodynamic (MHD) turbulence. We find good agreement for P{Ad) obtained in 
the two cases when the magnetic guide- field to fluctuations ratio Bo/brms is chosen accordingly. 
However, the scale invariance of P{AB/B) is broken in the MHD simulations with relatively 
limited intertial interval, which causes P{A6) to scale with measurement interval differently than 
in the solar wind. 

Subject headings: solar wind, magnetohydrodynamics (MHD), turbulence, interplanetary medium, plas- 



1. Introduction 

Since the beginning of the Space Age, the solar 
wind has been used as a natural laboratory for 
studying plasma turbulence. It is now widely ac- 
cepted that the s olar wind contains coherent mag- 
netic structures ([Veltri 1999 : Sorriso-Valvo et al 



1999l:lBruno et al.l2007ailLi et al.ll20f lllMiao eTal 



20111 ). Observations have shown that magneti- 
cally dominated structures are ubiquitous in the 
heliosphere and are advected by both the fast 



and slow wind (jBruno et al.ll2007bf ). There have 



been ma ny approaches to studying these struc 



tures (H l2008t ICreco et al.l [ioTot Iservidio et aD 



20111 ). hinging on the fact that they are typically 
associated with rapid spatial variation or rever- 
sal of the local magnetic field. In other words, 
structures are associated with discontinuities in 
the magnetic field. 



Although it is accepted that magnetic struc- 
tures exist in the solar wind, there remains a 
debate regarding the origin and nature of these 
structures. Two primary theories have been put 
forth. One model envisions that the structures are 
fiux tubes that originate in the solar corona and 
are passively advected by the solar wind. A sec- 
ond model contends that magnetohydrodynamic 
(MHD) turbulence dynamically produces inter- 
mittent structures as the solar wind fiows outward. 
Although these two models are not mutually ex- 
clusive, it is important to determine to what de- 
gree each mechanism contributes. 

An critical step toward answerin g this ques- 
tion was taken by iBorovskv ( 2008 ). He noted 
that magnetic discontinuities are implied by large 
fluctuations in the magnetic field direction across 
short measurement scales ( Bruno et al.l 12001 
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l2004l) . which provides a convenient and statis- 
tically robust way to quantify magnetic discon- 
tinuities. Borovsky therefore considered fluctu- 
ations in the magnetic field direction, given by 
rotations A6' = cos~^ (Si • B2/I-B1IIB2I) where 
Bi and B2 are the magnetic field vectors for 
two measurements separated by time At. Us- 
ing measurements from the ACE spacecraft, he 
studied the statistical properties of rotations 
by constructing a probability density function 
(pdf), P{A9). Borovsky discerned two popula- 
tions in the pdf. The first population consists 
of strong discontinuities at 30° < A9 < 170° 
with an exponentially decaying pdf proportional 
to exp(— A6'/24.4°). The second population con- 
sists of weak fiuctuations at 5° < A9 < 30° 
which can be fit by exp(— A 6'/9.4°). Subsequent 
anaylsis by iMiao et"aL ( 20 111 ) using data from the 
Ulysses spacecraft was in broad agreement with 
Borovsky's result. The interpretation put forth by 
these authors was that the first population rep- 
resents strong discontinuities across the walls of 
coronal flux tubes, while the second population 
represents turbulent fluctuations. 

However, it is known that MHD turbulence can 
spontaneo usly produce strong r nagnetic discon- 
tinuities (iGoldstein et all Il995l: iBoldvrevI l2006t 



Greco et al.l 120081 ). These discontinuities cor- 
respond to intermittent structures, which pri- 
marily take the form of current sheets for iner- 
tial range turbulence. It has also been demon- 
strated that many statistical properties of the so- 
lar wind fluctuations are consistent with MHD 
turbulence. These include, for instance, the 
energy spect ra of magnetic and velocity field 



fluctuations (jPodesta et al.l 120071 : iBoldvrev et al 



201 II ). the anisotropic s caling of structure func- 
tions (jChen et al.l 120111) . an d the waiting-time 
analysis of field increments ( Greco et al. 20091 



aiyi 

20101 ). The MHD description is known to break 



down at small scales, but the solar wind fluctua- 
tions are then thought to be consistent with ki- 
netic models of tu rbulence (jMedvedev et al.lll997 
Howes et al. 201ll ) . We have also shown in a sepa- 



rate paper that MHD turbulence is able to repro- 
duce the observed exponential tail of P{A9) for 
reasonab le strengths o f the magnetic guide-field 
jziidan kin et al. 2012). These facts support the 
picture where the population of discontinuities 
is governed by MHD turbulence, which may be 



co mplementary or alternative to inte rpretations 
bv lBorovsk^ (|2008h : [Miao et all (|201l[ ). 

In this Letter, we investigate the statistical 
properties of magnetic discontinuities in the so- 
lar wind. We show that P{A6) is a consequence 
of the magnetic field vector having a strong ten- 
dency to undergo pure rotations with lognormally 
distributed relative increments AB/B. This pic- 
ture is attractive because the magnetic field in the 
solar wind is known to exhi bit lognormal statis- 
tics to a good approximation ( Burlaga]|2001 ). Fur- 
thermore, we find that the pdf 's of relative incre- 
ments, P{AB / B), for different separations At can 
be rescaled to match, that is, the distribution of 
{AB/B) ■ (Ai/Ato)^" with appropriately chosen 
AIq and a is a universal lognormal distribution. 
This strong scaling property implies the so-called 
monofractality, or absense of intermittency in the 
fluctuations of AB/B. 

We then directly compare the statistical prop- 
erties of rotations in the solar wind with numer- 
ical simulations of MHD turbulence. We flnd a 
good agreement in P{A6) for the two cases for 
reasonable ratios of simulation guide-field to fluc- 
tuations. We also flnd that the scaling of the dis- 
tribution function of AB/B is not well preserved 
in the simulations and that the magnetic fleld vec- 
tor does not have a strong tendency to undergo 
pure rotations. These differences may be the con- 
sequence of a much smaller inertial range in nu- 
merical simulations than in the solar wind. This 
causes the pdfs in the two cases to scale differently 
with measurement interval, and hence separations 
in the two cases cannot be directly compared. 

2. Analysis of solar virind 

First, we consider the properties of magnetic 
discontinuities in the solar wind. We analyze a 
time-series of spacecraft magnetic field measure- 
ments, B{t), taken by the Wind spacecraft over 
the years 2005-2010. Each vector is averaged over 
a time of Atmin = 3 seconds, but the results 
are independent of the averaging scale as long 
as At > Atmin- The analysis does not exclude 
events such as CMEs and CIRs. It also does not 
differentiate between slow and fast winds, which 
are known to have different typical properties and 
physical processes. A similar statistical study ap- 
plied separately to the regions of fast or slow wind 
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may in principle produce different results. 

To quantify magnetic discontinuities, we con- 
sider fluctuations in the magnetic field direction, 
given by the rotation A^. The rotation /S.0 be- 
tween two magnetic field vectors, Bi and B2, sep- 
arated by time interval At, is 



cos 



(61 • b2 



(1) 



where b = B/B is a unit vector along the mag- 
netic field. A large rotation implies a magnetic dis- 
continuity, while a small rotation represents weak 
fluctuations. The statistical properties of rota- 
tions then provide a window into understanding 
the nature of discontinuities and intermittency. 
Similar measurement s of magnetic field rotations 
have been applied bv IPerri et al. ( 2009[ ) to study 
anisotropy. 

In Fig.[TJ the pdf of rotations, P{A9), is shown 
for At e {1, 16, 256, 4096} • 3 seconds. The pdf has 
a strong dependence on At. The tail {A9 > 60°) 
can be approximately fit by an exponential func- 
tion, while the bulk of the pdf {Ad < 60°) varies 
drastically. This reveals that the resolution used 
to measure fluctuations is an important param- 
eter. This is reasonable because the correlation 
between measurements is increased when At is de- 
creased, which in turn increases the likelihood of 
small A^. 

We now show that P{A9) can be described by 
a simple model. We define the relative increment 
in magnetic field as AB/B = \B2 - -Bi|/|-Bi|. If 
the magnetic field vector undergoes a pure rota- 
tion between the measurements, that is, does not 
change amplitude, then AB/B depends only on 
the rotation angle A9, satisfying 



AB/B = I62 - bil = 2sin(A6l/2), 



(2) 



where fluctuations must lie within < AB /B < 2 
to be meaningful. We find that the magnetic field 
vector in the solar wind indeed strongly satisfies 
Eq. [21 as shown in Fig. [51 Therefore, between two 
closely spaced measurements {At < 10^ seconds), 
the magnetic field tends to mostly rotate instead 
of change in strength, and the fiuctuations in the 
rotation angle can be studied in terms of the rela- 
tive increments AB/B. 

We find that P{AB / B) in the solar wind is well 
fit by a lognormal distribution for < AB/B < 2. 
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Fig. 1. — The pdf of rotations in the magnetic field 
direction, P{A9), for the solar wind with measure- 
ment intervals of At £ {1, 16, 256, 4096} -3 seconds 
(measured by the WIND spacecraft). While part 
of the tail {AO > 60°) may be fit by an exponen- 
tial function (straight line in the plot), the bulk of 
the pdf {A9 < 60°) varies drastically with At. 
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Fig. 2. — A scatterplot of 2sin(A6'/2) versus 
AB/B, which are equal if the magnetic field 
vector undergoes a pure rotation between mea- 
surements. Also shown, in the inset, is the 
pdf of the normalized difference x = \AB/B — 
2sm{A9/2)\/{AB/B), which is strongly peaked 
at zero. These plots show that the magnetic field 
in the solar wind tends to rotate between the mea- 
surements rather than change in strength. 
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The lognormal distribution is given by 



/(^) 



1 



xay2'K 



cxp 



(3) 



where and a are the location parameter and 
scale parameter, respectively. The fit is not sur- 
prising since other measurements of magnetic 
fields in the solar wind also have appr oximately 
lognormal statistics, e.g. iBurlaga ( 2001 ). The log- 
normal fits to P{AB/B) can all be chosen with 
(7 = 1, which implies that P{AB/B) for any At 
can be rescaled into the universal lognormal func- 
tion according to the formula 



At 
Ak> 



F 



( AB f At_ 



(4) 



where F{x) is the universal lognormal probability 
density function, 



exp ( — - lo. 



(5) 



and the best fit is given by Ato ~ 6.6 x 10^ sec- 
onds and a w 0.46. This property is illustrated 
in Fig. m where the pdfs of the rescaled vari- 
able X = {AB/B) ■ (At/Aio)"" are in remarkable 
agreement for 3 < At < 1024 • 3 seconds. We 
found that the measured curves follow the scale- 
invariant lognormal form almost until the cutoff 
that occurs at AB/B ~ 2. For different values 
of At such a cufoff corresponds to different values 
of X. Fig. [3] shows the range of x, where all the 
measured curves follow the lognormal form F{x). 

This scaling property of the probability density 
function implies that the field AB/B has a scale- 
invariant and nonintermittent (monofractal) dis- 
tribution. In addition, it demonstrates that there 
exists a certain scale Ato-, which can probably be 
identified with the outer scale of turbulence. 

We therefore model the rotations Ad as fol- 
lows. Assume that P{AB/B) is lognormal and 
scales with At as described by and also as- 
sume that the vector undergoes a pure rotation so 
^ is apphcable. Then P{Ae) is obtained. This 
model gives excellent agreement with the pdf of 
rotations in the solar wind, as shown in Fig. |4] for 
At e {1, 16, 256, 4096} • 3 seconds. It is worth not- 
ing that the bulk of P{A9) is also approximately 
lognormal, in agreement with previous conclusions 
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Fig. 3.— The quantity x = {AB/B) ■ (A</Ato)~" 
has the pdf F{x) that is independent of At, where 
Ato ~ 6.6 X 10^ seconds and a « 0.46. The curves 
shown for At G {1, 4, 16, 64, 256, 1024} • 3 seconds 
are all in remarkable agreement with the lognor- 
mal distribution of = and a ~ I, indicating 
that AB/B is scale- invariant over the given values 
of X. For small At, the agreement also extends far 
beyond the values of x shown, almost up to the 
cutoff at AB/B = 2. 
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by IVasQuez et all (|2007l ) and iBruno et all (|2004[ ). 

We have however shown that modehng Ai?/B, 
rather than A0, from a lognormal distribution bet- 
ter fits the data. The weak fluctuations and strong 
discontinuities characterized by /S.9 may then be 
regarded as a single population described by log- 
normal fluctuations in the magnetic field and a 
tendency of the magnetic field vector to rotate. 

3. Analysis of MHD turbulence 

We now discuss the statistical properties of ro- 
tations in MHD turbulence, and compare them to 
the solar wind case. We use data from direct nu- 
merical simulations of 3D incompressible strong 
MHD turbulence. The MHD equations are 



dtv + [v ■ \7)v 
dtB 

V • V 

V B 



-Vp + (V X B) X J5 + i/V^f + fi, 

V X (i; X B)+?7V2jB + f2, 

0, 

0, (6) 



where v{x,t) is the plasma velocity, B{x,t) = 
Bq -I- b is the magnetic field that has both uni- 
form (Bo) and fluctuating (b) components, p is 
the pressure, and f(a;,i) is the external forcing. 
We take the viscosity v and resistivity "q to be 
equal. 

The equations are solved on a triply periodic 
domain using standard pscudospcctral methods. 
The time-advancement of the diffusive terms is 
carried out exactly using the integrating factor 
method, while the remaining terms are treated 
using a third-order Runge-Kutta scheme. For a 
detailed d escription of the nume rical method, see. 



e.g., Ref. ( Cattaneo et al. The turbulence 



is driven at the largest scales by applying random 
forces fi and £2 in Fourier space at wave-numbers 
2tt/L < k± < 2{2tt/L), fc|| = 27r/L. The corre- 
lation between the forces is chosen as to mimic 
the driving by independent co unter-propagating 
shear - Alfven modes, see, e. g., (jPerez fc BoldvrevI 
2OIOI : Boldvrev et al. 2011 ). however turbulence 
may also be driven by other choices of large- 
scale forcin g; this does not aff ect the inertial in 



terval, see ( Mason et al. 2008[) . The forces have 



no component along z and are solenoidal in the 
a;y-plane. All of the Fourier coefficients outside 
the above range of wave-numbers are zero and in- 
side that range are Gaussian random numbers that 



are refreshed on average every Q.lL/[2-KVrms) time 
units (that is, force is updated approximately 10 
times per large-scale turnover time) with ampli- 
tudes chosen so that Vrms ^ 1- We conducted a 
number of MHD simulations with different ratios 
of Bo/brrns- The analysis is performed on several 
statistically independent snapshots corresponding 
to a steady-state, with resolution of 1024'^ and 
Reynolds number Re ~ 3000. 

For a fixed time snapshot, we measure the ro- 
tation of the magnetic field between points sepa- 
rated by spatial distance Ax. Therefore Bi and 
B2 in Eq.[l]are taken from a magnetic field profile 
B{x,y,z). We consider separations taken in the 
plane perpendicular to the guide-field, although 
the results are similar for separations along the 
guide-field. Using spatial separations in MHD in- 
stead of time separations is a significant difference 
from the solar wind analysis. However, we ex- 
pect temporal measurements in the solar wind to 
be equivalent to spatial ones because of the Taylor 
hypothesis. At « Ax/Vswi where Vsw is the (ap- 
proximately constant) solar wind velocity. There- 
fore, we should observe At oc Ax if the solar wind 
observations are consistent with MHD turbulence. 

We find that P{A6) has a strong dependence 
on the ratio of magnetic guide-field to root mean 
square (rms) fiuctuations, Bo/b„ns- Specifi- 
cally, the characteristic angle for the exponen- 
tial tail decreases with increasing Bo/hrms^ i-e. 
the slope on the log-h near axis becomes steeper 



(jZhdankin et al.l l20f 21 ). We find that choosing 



Bo/brms ~ 0.32 gives a good agreement with the 
solar wind measurements considered here, which 
is reasonable since Bq &rms in the solar wind. 
Therefore, we use MHD data with Bo/brms ~ 0.32 
in the remaining analysis. We also find that 
P{A6) does not change significantly in simula- 
tions of lower Reynolds number {Re « 2200) 
and lower resolution (512"^), indicating that even 
higher Reynolds numbers may need to be achieved 
for better agreement with the observational data. 

A direct comparison of P{Ad) for MHD tur- 
bulence and for the solar wind is shown in Fig. [5] 
for several measurement intervals. To get the best 
agreement between the two cases, we have cho- 
sen intervals of Ax/L € {4, 16, 64}/1024 for MHD 
and At € {5, 64, 830} • 9 seconds for the solar wind. 
This requires approximately that At oc Ax^'^^ in 
order to get the solar wind pdfs to match the 
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Fig. 4. — Comparison of the pdf for rotations mea- 
sured in the solar wind (in blue) and the pdf de- 
rived from lognormal P{AB/B) with the assump- 
tion that vectors undergo pure rotations (in red). 
The separations shown are At e {1, 16, 256, 4096}- 
3 seconds. 



MHD pdfs, which seemingly violates the Taylor 
hypothesis. We find also that the scale invariance 
of P{AB/B) as described by (jlj is broken in our 
MHD simulations: although P{AB/B) is well fit 
by a lognormal distribution for all At, the fits re- 
quire different values of a and therefore cannot be 
rescaled into one another. We also find that the 
magnetic fluctuations in the simulations are not 
dominated by rotations. 

These differences are likely caused by a much 
shorter inertial range (less than one decade) that 
spoils the scale invariance of P{AB/B) as de- 
scribed by (|4]), which in turn modifies the scal- 
ing of P{Ad). One, in princple, should not expect 
to see such a strong scaling property of AB/B in 
a rather modest simulation. A related difference 
is that the MHD measurements are near the dis- 
sipation range (Ax < 16), while the solar wind 
measurements are well within the intcrtial range. 
Therefore the simulations must have a larger sep- 
aration of scales to more accurately capture the 
physics that describes the scaling. Nevertheless, 
we see that for any fixed separation, MHD is able 
to produce similar pdfs as in the solar wind, which 
shows that discontinuities in the solar wind are ac- 
curately reproduced by simulations. 
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Fig. 5. — A comparison of the pdfs of A9 for MHD 
turbulence (in red) with Ax/L € {4, 16,64}/1024 
and for the solar wind (in blue) with At G 
{5, 64, 830} • 9 seconds, consistent with the relation 
At ^ (Aa;)^-*'^. There is generally good agreement 
in the pdfs for both cases. 



4. Conclusions 

We found that the pdf of rotations in the mag- 
netic field of the solar wind has a shape that is in 
good qualitative agreement with strong MHD tur- 
bulence (Fig. [5]). The numerical simulations are 
consistent with the observations if the magnetic 
guide- field to rms fluctuations ratio is Bo/brms ~ 
0.32, which is reasonable for the solar wind. The 
approximately exponential tail of P(A 9) is con- 
sisten t with solar wind a nalysi s done bv iBorovskv 
(|2008l ) and iMiao et al.l (|201ll) . but our interpre- 
tation of the result is different. Our results sug- 
gest that the pdf is closely associated with log- 
normal fluctuations of the magnetic field. Since 
such fluctuations are present in MHD simulations, 
this suggests that MHD turbulence can, to a large 
extent, describe the magnetic discontinuities ob- 
served in the solar wind. This supports a pic- 
ture in which the discontinuities are associated 
with intermittent structures that arise from the 



MHD energy cascad e (iSorriso-Valvo et al.l 12007 ; 
MacBride et al.ll2008l ). This can be alternative or 
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complementary to the coronal flux tube model of 
discontinuities. 

What remains to be better understood is the 
difference in scaling of Ax and At required to ob- 
tain agreement between the pdfs for the two cases. 
This result is surprising, and wc cannot provide 
an exhaustive explanation with our current level 
of understanding and numerical simulations. Our 
results may suggest that the different scaling in 
MHD simulations arises due to a very limited in- 
ertial range relative to the solar wind. The fact 
that P{A9) are in such good agreement despite 
these differences suggests that there could be an 
underlying universal description for them. 
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